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We have developed a model to describe the interactions of neutrinos with nucleons and nuclei, 
focusing on the region of the quasielastic and A(1232) peaks. We describe neutrino-nucleon collisions 
with a fully relativistic formalism which incorporates state-of-the-art parametrizations of the form 
factors for both the nucleon and the N — A transition. The model has then been extended to finite 
nuclei, taking into account nuclear effects such as Fermi motion, Pauli blocking (both within the local 
density approximation), nuclear binding and final-state interactions. The in-medium modification 
of the A resonance due to Pauli blocking and collisional broadening have also been included. Final- 
state interactions are implemented by means of the Boltzmann-Uehling-Uhlenbeck (BUU) coupled- 
\^**) ' channel transport model. Results for charged current inclusive cross sections and exclusive channels 

as pion production and nucleon knockout are presented and discussed. 
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i. introduction 
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Neutrino experiments now provide conclusive evidence that neutrino oscillations exist an( i; therefore, neutrinos 

£N| are not massless. An extensive experimental program aiming at the precise determination of v masses, mixing 

angles and CP parameters is being developed in several labs around the world. The success of these oscillation 
experiments depends critically on a good knowledge of neutrino-nucleus interactions. In this way it is possible to 
minimize the systematic uncertainties in neutrino fluxes, backgrounds and detector responses. 

The interest in neutrinos goes beyond the study of their intrinsic properties and extends to a variety of topics 
in astro-, nuclear and hadronic physics. In particular, they are a valuable tool for nuclear and hadronic physics. 
\Q , The availability of a high intensity v beam at Fermilab offers a unique opportunity to gain new information on the 
structure of the nucleon and baryonic resonances. Experiments such as MINERVA |(| and FINeSSE .7] shall address 
relevant problems like the extraction of the nucleon and TV — A axial form factors, or the measurement of the strange 
spin of the nucleon However, those experiments will be performed mainly on nuclear targets since they provide 
relatively large cross sections. 

Particles produced in neutrino interactions can reinteract before leaving the nucleus and can be absorbed, change 
their kinematics or even charge before being detected. For example, even if an initial interaction might be quasielastic 
scattering the observed final particles can include pions and multiple nucleons. Thus, understanding nuclear effects 
is essential for the interpretation of the data and represents both a challenge and an opportunity |jj ■ 

There is a rather extensive literature on the theoretical study of neutrino- nucleus collisions. Nuclear shell model [To| 
and RPA calculations [Til IT^] address the region of low energy neutrinos. Those reactions are important for the 
d ' synthesis of heavy elements in stars 0, 0] . 

At intermediate energies, quasielastic collisions have been investigated using a relativistic treatment of the matrix 
element in Fermi gas models |0, 0, that take into account Fermi motion, Pauli blocking and binding energies. 
The importance of strong renormalization effects on the weak transition strengths in the nuclear medium has also 
been emphasized. These nuclear correlations are often taken into account as RPA resummation of particle-hole and 
A-hole states 0, 0, |2(j , but realistic nuclear spectral functions have been used as well [2l| • 

Nucleon knockout has attracted considerable attention, mainly in connection with the possibility of extracting the 
axial strange content of the nucleon [22l 123^ . The most straightforward approach to this problem is the plane wave 
impulse approximation. It neglects all interactions between the outgoing nucleons and the nucleus p4j . but it has 
been claimed that the distortions in the observables due to rescattering largely cancel for the ratio of proton to 
neutron yields |25|. This situation can be improved by means of a distorted wave impulse approximations, where the 
final nucleon scattering state is computed using optical potentials |26l |27| . For energies above 1 GeV the Glauber 
model, which is based on the eikonal and frozen spectator approximation, is a better alternative [2(|. However, in 
those approaches the nucleons that go into unobserved states as a result of an interaction are lost, while, in fact, 
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they are ejected with a different energy, angle and maybe charge. Monte Carlo methods permit to take into account 
interactions of nucleons leading to energy losses, charge exchange and multiple nucleon emission psj . 

Theoretical studies of neutrino-induced pion production in nuclei are considerably more scarce. The main pion 
production mechanism is via the excitation and subsequent decay of the A resonance p9l l30l l3l| . The contribution 
of higher order resonances at v energies < 2 GeV has been found to be small j3^,|3i|. Singh et al. |2^| pointed out 
that the in-medium modification of the A width causes a reduction of the v A — ► A X cross section. The final-state 
interactions of the produced pions is an essential ingredient of any realistic model of this reaction. However, the 
available calculations deal with it in an oversimplified way. Singh et al. only estimate the pion flux reduction in 
an eikonal approximation, neglecting elastic and charge changing scattering. In Refs. H3,|3j those processes are 
considered in a random walk approximation where the pion energy stays constant during their propagation out of 
the nucleus. The nucleons are frozen within the nucleus (no Fermi motion) and Pauli blocking appears only as a 
global factor common to all collisions. However, this approach does not take into account the knowledge gathered in 
extensive research of pion production in nuclei using different probes in the past years |35| . 

Finally, at higher energies, deep inelastic scattering becomes progressively important l3fl. The role of nuclear 
corrections in this case is currently debated in connection with the recent NuTeV findings |37j . The impact of deep 
inelastic scattering in the energy region under study here is expected to be minor. Therefore, we do not take it into 
consideration. 

Quasielastic scatterin g, p ion production and deep inelastic scattering are put together in event generators (e. g. NU- 
ANCE [38|. NEUGEN |39| and NEUT ^1) used by the neutrino experiments in their simulations. For resonance 
production, these generators take the model of Rein and Sehgal |4l| based on old quark model calculations for the 
form factors. They differ substantially in the implementation of in-medium effects and final-state interactions. For 
example, NEUT and NUANCE apply a relativistic Fermi gas model with a global Fermi momentum and a constant 
nuclear binding energy. The in-medium modification of the A width is taken into account only in the NEUT model. 
The final-state interactions are treated by using a cascade model with rescattering, charge exchange and absorption 
of pions. 

In this work we study neutrino interactions with nuclei at intermediate energies where we have included the two 
most important contributions to the vN collision at neutrino energies up to about 1.5 GeV, namely quasielastic 
scattering and A excitation. As in photonuclear reactions we assume that the neutrino interacts with a single bound 
nucleon in the nucleus. Thus, the first essential step is a reliable description of the vN collision. For that we use 
a fully relativistic formalism with state-of-the-art parametrizations of the form factors. Then, the elementary cross 
section has to be modified in order to account for nuclear effects. We describe the nucleus by a local Fermi gas, 
leading to density dependent Fermi momenta. We further account for the binding of the nucleons in a density and 
momentum dependent mean-field. Pauli blocking of the final state nucleons is also considered. 

Once produced, the particles experience final-state interactions when propagating through the nuclear medium. A 
realistic treatment of the final-state interactions can be achieved in the framework of a coupled-channel transport 
theory — the Giessen BUU model — which allows the investigation of exclusive channels. Our results show that the 
coupled-channel effects influence in particular the low energy part of the hadron spectra by side-feeding — this is an 
important feature of our model that cannot be achieved with simple absorption models, as e. g. Glauber approaches. 
Furthermore, we emphasize that our transport model has been extensively and successfully tested against experimental 
data for 7^4, eA, ttA and also heavy ion collisions. Thus, the same theory is able to describe many different nuclear 
reactions with the same set of parameters and physics assumptions. Since there are no free nuclear parameters for 
neutrino-nucleus reactions, our theory should be able to give reliable predictions for vA scattering relevant for future 
neutrino experiments. 

This article is organized as follows. In Sec. [n] we first present our model for neutrino-nucleon scattering. Then we 
discuss the in-medium modifications needed in the case of scattering with a bound nucleon. We give a brief overview 
of the BUU transport model that is used for the calculation of exclusive channels. Our results for inclusive cross 
sections, pion production and nucleon knockout are presented in Sec. IIIII Our model is general since it includes all 
v flavor for both neutrino and antineutrino, as well as charged and neutral current in any arbitrary nuclei. Here , w e 
focus on charged current Vy, scattering off 56 Fe. First results of our model were already presented in Ref. p^l43l l44 | . 
We close with a summary and conclusions in Sec. IIVI 

II. THEORY 

In the impulse approximation, the neutrino-nucleus interaction can be treated as a two-step process. The first 
step is the interaction of the neutrino with a single bound nucleon in the target nucleus. Then, as second step, we 
propagate the produced particles through the nucleus using the BUU transport model including a full coupled-channel 
treatment of the final-state interactions. Below we describe those ingredients in some detail. 
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FIG. 1: Total neutrino cross sections for v^N — > yT X for an isoscalar target as a function of the neutrino energy decomposed 
in the contributions of quasielastic scattering (QE), A excitation (RES) and deep inelastic scattering (DIS). The DIS curve is 
taken from Ref. |3q|: the QE and RES curves are calculated with the formalism presented in Sec. Ill Al averaging over protons 
and neutrons: ct^jv = {(Tu^p + o"i/ (1 n)/2. 



A. Elementary neutrino-nucleon reaction 



The cross section for the elementary interaction of a neutrino with a nucleon, 

v{k)+N{p)^l-{k')+X(p% (1) 

with the momenta k a = ^E v , fc V k' a = (Ei, k'^j , p a — (E, p) and p' a = [E' , p'), is given by 



f d 4> 1 9i; 1 ±s(p' 2 -M'' 
J v 64tt 2 \k-p\ E v V 



\M\ 2 (2) 



&Q 2 m 

with Q 2 = — (k — k') 2 and the azimuthal angle <j> £ (0,27t). M' denotes the (invariant) mass of the outgoing baryon. 
The matrix element squared \M\ 2 for charged current (CC) interactions, summed and averaged over spins, reads 

\M\ 2 = ^f^L aP H^, (3) 

where the leptonic tensor, L a f3, is given by 

L aP = Tr[(^+m !/ )7 a (l-7 5 )(^+m07^(l-75)] 

= 8 [k' a k (j + k a k' p - g a(3 k ■ k' + ie^p^k" 7 ] . (4) 

mi is the mass of the outgoing lepton. It has been shown, that in particular at low momentum transfer, the lepton 
mass becomes important l4rij. Thus, we do not neglect it like in other works (e. g. Ref. The neutrino mass 

m u , however, is set to zero. For the corresponding antineutrino reaction, 

D(k)+N(p)^l+(k')+X(p'), (5) 

the antisymmetric term proportional to the fully antisymmetric tensor e a p p(7 gets a minus sign. 

The hadronic tensor H a P depends on the specific reaction. In the intermediate energy region, i. e. up to neutrino 
energies of about 1.5 GeV, the reaction is dominated by two processes, namely quasielastic scattering (QE) and A 
production (cf. Fig. Restricting to CC interactions, we consider in this work 

vn — * l~p, (6) 

(note that the charge changing QE reaction is not possible on protons) and 

wnrA+, vp~>rA ++ . (7) 
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Both contributions to the vN cross section are described within a fully relativistic formalism which is presented below. 

We want to note that, even though quasielastic reactions and A production are the most important processes in the 
energy regime under investigation, for a full description one has to include higher-mass resonances, the non-resonant 
background and deep inelastic scattering (DIS) as well as more exotic channels such as strangeness production. 



1. Quasielastic scattering 

Since quantum chromodynamics (QCD) does not allow a direct calculation of the hadronic vertex, we need to 
construct an explicit expression for the hadronic current of the QE reaction based on general assumptions. Following 
the arguments of Ref. 47] we construct the most general form from the four- vectors at our disposal: p a , p' a and 
Qa = Pa ~ Pa- Gordon identities (cf. e. g. Appendix A. 2 in Ref. |48|) limit the number of terms and we obtain 



J% E = (p\JS E (0)\n) 
= u(p')A a u{p) 



with 



f ) rt 



2M N 



a a fsq^F^ + 7 Q 75 -Fa 



M N 



75 F F 



(8) 



(9) 



with Mn being the nucleon mass. F^ 2 are the vector form factors; Fa and Fp are the axial and pseudoscalar form 
factor, respectively. The term {Jiq a )/q 2 appears when the masses of in- (M) and outgoing nucleons (M') are not 
equal. It ensures that the vector part of the current is conserved even for nonequal masses. This is the case in the 
nuclear medium where a momentum dependent mean-field can cause such a difference in masses (cf. Sec. Ill B|) . This 
extra term vanishes by applying Gordon identities in the limit of M — M' , and Eq. @ reduces to the "standard 
expression" (see e. g. Ref. |49|) 



-4, 



2M N 

The hadronic tensor iJ Q/3 for QE reactions is given by 

1 



Tra/3 
QE 



-Tr 



o a(i q Si F 2 v + "f a "f b F A + Ty-75-Fp- 
M N 



{i>+M)A a {i>'+M')A^ 



(10) 



(11) 



with 



A n 



(12) 



and — in the free nucleon case — with M — M' — Mm- 

The vector form factors F^ 2 can be related to electron scattering form factors. The conserved vector current 
hypothesis (CVC) implies that the vector part of the current in Eq. and the electromagnetic current are components 
of the same isospin multiplet of conserved currents, and that therefore their form factors are related. We obtain 



fY{Q 2 ) 



{G p E - Gl) + 



Q 2 



4M N 



Q 2 



AM N 



(13) 
(14) 



Gm(Q 2 ) and Ge(Q 2 ) are the magnetic and the electric form factors of the nucleon respectively, for which we take the 
BBA-2003 parametrization [5(j. It uses recent electron scattering data from JLab to account for deviations from the 
dipole Q 2 -dependence. By assuming pion pole dominance, we can use the partially conserved axial current hypothesis 
(PCAC) to relate F A and F P , 



F P (Q 2 



2M% 



-F A {Q 2 



(15) 



The axial form factor is given by the standard dipole form 

Fa{Q 2 ) = 9 A 



1 



M 2 A 



(16) 
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FIG. 2: Total QE cross section as a function of the neutrino energy (solid line). The datapoints are taken from Refs. |5l| (•), 
52] (■), |H^| (A). The dashed lines are the results of a 5% change in Ma in Eq. 1161 . 



the axial vector constant qa — —1.267 is obtained from (3 decay. The axial mass Ma has to be measured in neutrino 
experiments. We use Ma = 1-00 GeV which was extracted from data by the analysis of Budd et al. |Hcj . 

With these ingredients we obtain the QE cross section on the neutron, which is presented in Fig. [21 in comparison 
with the available data on H2 and D2 targets. The main source of uncertainty comes from the axial form factor 
whose Q 2 dependence has to be extracted from v scattering data. Experimental analyses (see table 3 of Ref. [5(| and 
references therein) give Ma with an accuracy of about 5%. The impact of this error on the cross section is shown by 
the dashed lines. It is clear that more precise measurements are required to constrain more this parameter. 



2. A production 



For a theoretical treatment of the neutrino-induced charged current A production we use a fully relativistic formal- 
ism (cf. Refs. j2}||22,E10.IE3)- ^he hadronic current for the reaction in Eq. JJJ is given by 



with 



C v 



J A = <A+|J Q A (0)|n) 



j^ridam ■ p' - qpp'a) + j^ridapq ■ p - qpp a 



75 



M N 



qA qA 

+ T7-(fa/3^ ~ qpla) + -Tjr(g a pq ■ P' ~ qpP' a ) + C£g a f) + -r^-qfSqa, 
1V1 N IV1 N 



(17) 



(18) 



where ip^{p') is the Rarita-Schwinger spinor for the A, and u(p) is the Dirac spinor for the nucleon; q — k — k! is the 
momentum transfer. This yields the following hadronic tensor 



H 



a/3 



-Tr 



# + M)B a PA pa B al3 



with 



B a p = 70^70 

and, for free nucleons, M = Mm- The spin 3/2 projection operator is given by 



A, 



9pa 



2 PpP* 

3 p' 2 



1 P'plo ~ Palp 1 



3 7p7a 



(19) 



(20) 



(21) 
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By isospin relations, we obtain for A ++ production: 

(A++\J a (0)\p) = \/3<A+|J Q (0)|n>. (22) 

Basically, two approaches are considered in the literature for the N—A transition form factors C^' A with i = 3, . . . , 6. 
They can be phenomenological with parameters extracted from neutrino and electron scattering data or calculated 
using quark models. Early attempts for the latter are reviewed in Ref. [54| (see also references therein); more recent 
ones are summarized in Ref. |5(| . In their model of resonance production, Rein and Sehgal |4l| adopted the quark 
model of Feynman, Kislinger and Ravndal 153 ■ A more recent calculation was done by Liu et al. |5a ] who applied the 
Isgur-Karl quark model and by Sato et al. [58j who developed a dynamical model including pion cloud effects. 

As in QE scattering, we choose the first approach and use phenomenological form factors. The vector form factors 
CY can be related to the ones obtained in electron scattering by applying CVC. This, together with the assumption 
of M\ + dominance of the electroproduction amplitude, yields [32J 

C^(Q 2 ) = 0, C 5 v (Q 2 )=0 and C\ {Q 2 ) = — ^=C^(Q 2 )- (23) 



This leaves only one independent vector form factor, C% , which can be parametrized in various ways to describe 
pion electroproduction data. The data show that the Q 2 -dependence is steeper than that of a dipole. Thus we use a 
modified dipole form and adopt a parametrization used in Ref. [33] ] 

W ) = C ?(0)(l + $)"(l + J ^)", P4) 

with CY(0) = 1.95 and M v = 0.84 GeV. 

In the axial sector, we apply similar techniques as in the QE case. Pion pole dominance yields for C§ [54| 

C ^>= X^V ' ( ° 2) (25) 

with being the pion mass . gAN-w is the A ++ — > prr + coupling constant and / T , the pion decay constant. Ft T (Q 2 ) 
is the vertex form factor with F^(Q 2 — m 2 ) = 1. This relation together with the assumption of PC AC connects C§ 
and C§ in the axial current 

C £(Q 2 ) = Ci{Q 2 )-^. (26) 

In the limit Q 2 = and with the assumption that F^(Q 2 ) is a slowly variating function with F 7T (Q 2 — m 2 ) i=s F 7T (Q 2 = 
0) = 1, we obtain the off-diagonal Goldberger-Treiman relation 

C 5 A (0) = 9A "*f* ~ 1.2. (27) 

This coupling was extracted from the BNL data in Ref. ^5| and found to be consistent with the PCAC prediction. 

Since there are no other theoretical constraints for C^(Q 2 ), Cf(Q 2 ) and C^(Q 2 )/C^(0), they have to be fitted to 
neutrino scattering data. The available information comes mainly from two bubble chamber experiments, the 12-foot 
bubble chamber at Argonne ( ANL) [5^. l60| and the 7-foot bubble chamber at Brookhaven (BNL) [6l| ■ The fits 



adopted the Adler model |62| where 



C£{Q 2 ) = -91^1 and C 3 A (Q 2 ) = 0. (28) 



4 

For we have taken, as in Ref. [3^, again a modified dipole 

C ^)= Ct »(l + ^)"(l + ^)", (29) 

with Ma = 1-05 GeV. We conclude that the neutrino- induced A excitation is fully described by two independent 
form factors, C\ and Cg . 
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The finite width of the A resonance is accounted for in the cross section of Eq. J5J by replacing 

S (p /2 - M 2 A ) -> - X -Lm I - — - . __ ] , (30) 



where we used M' = M A with M A = 1.232 GeV. 

For the free decay width T we use an energy dependent P-wave width as required by angular momentum conserva- 
tion. It can be parametrized by using Blatt-Weisskopf functions 

r - r °W (31) 

with T = 0.118 GeV, W = and 

m) ~ w i + { qcM {w)Ry {i2) 

For the interaction length R we take R — 1 fm. <7cm is the pion momentum in the rest frame of the resonance, for 
which we obtain 

qcM{W) — • 1^3) 

With these ingredients, we obtain the cross section for CC A ++ and A + production on proton and neutron, 
respectively. Once produced, the A decays into irN pairs. For the amplitudes A of pion production the following 
isospin decomposition applies 

A{v lP ^r P n + ) = A 3l (34) 
1 2 /2 

A{vin -> rnn + ) = -A 3 + —A u (35) 

A(nn^r P ir°) = -^-A 3 + -A!, (36) 

with A3 being the amplitude for the isospin 3/2 state of the irN system, predominantly A, and Ai, the amplitude 
for the isospin 1/2 state which is not considered here. Those three channels are plotted for in Fig. [3] together 
with pion production data from ANL and BNL on H2 and D2 targets. Since 7r + p production is dominated by A 
excitation, this channel serves as a quality check of our model, and a good agreement with the data is reached. In 
the case of 7r + n and 7r°p our results lie systematically below the data. This is mainly due to the non-negligible 
contribution from the heavier isospin 1/2 resonances like JV(1440), iV(1520) and iV(1535), as can be seen in Ref. 
where those states are considered. We achieve a good agreement with the ANL data when higher invariant masses 
are excluded (W < 1.4 GeV — open squares and open circles) |5^,|63- For the BNL experiment, data with cuts are 
not available, however, their analysis indicates, that at energies of 0.5 - 6 GeV approximately half of the 7r + n and 7r°p 
events correspond to invariant masses W > 1.4 GeV (cf. Fig. 6 of Ref. [(34j|). The underestimate in the (small) 7r + n 
channel translates into an error of about 10% in the isospin-averaged total 7r + production cross section at E u = 1 
GeV. This is well within any experimental uncertainties. 

As in the QE case, the main source of error for A production is contained in the axial form factors and, in particular, 
in which gives the dominant contribution to the cross section at low Q 2 . We have investigated the sensitivity of the 
cross section to the value of C§ at Q 2 — assuming an error of 5% as given in Ref. 45] (dotted lines). Furthermore, 
we show the impact of the uncertainties in the Q 2 dependencies of the form factors by assuming a 7% error in the 
determination of Ma as reported in the ANL and BNL analyses |(3(], EHI (dashed lines) . 

B. In-medium modifications of the uN reaction 

If the nucleon is not free but bound in the nucleus we must account for nuclear effects such as Fermi motion of the 
initial nucleons and Pauli blocking of the final ones. In local density approximation (LDA), the local Fermi momenta 
of the nucleons are given by 



3 



1/3 



?HO=(fV(r)] (37) 
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FIG. 3: Total pion production cross sections through A excitation as a function of the neutrino energy (solid lines) compared 
to the pion production data of Refs. |59| (• without and o with mass cut at W < 1.4 GeV), [g3 (B without and □ with mass 
cut at W < 1.4 GeV), |64|l (x), |6^| (A). Dashed and dotted lines reflect the uncertainties of C5 (0) and Ma, respectively. See 
text for the discussion of the discrepancies between the calculation and the data in the lower panels. Note also the different 
scales. 



For the density distribution we take a Woods-Saxon parametrization 

P{r) =Po (l + exp^— ^ , (38) 

with parameters extracted from Hartree-Fock calculations 66] . The corresponding values for 56 Fe are ro = 4.22 fm, 
a = 0.477 fm and po — 0.158 fm~ 3 . In neutrino-induced reactions, the phase-space density / is given by 

f(?,p)~e(p F (r)-\p\). (39) 

We can then approximate the probability that a final state of a nucleon is not Pauli blocked by 

P Pau u = l-e(p F (r)-\p\). (40) 

Furthermore, we take into account that the nucleons are bound in a density and momentum dependent scalar 
potential 

U = U(r,p), (41) 

which can be related to an effective mean-field potential as outlined in the following. The general expression for a 
relativistic one-particle Hamiltonian for a particle of mass m is given as 



H = 



sj{m + Sf + (p - U v ) 2 + K, (42) 
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where S is a scalar potential and (U v , Uy) a vector potential. The spatial components of Uy vanish in the rest frame 
of the nucleus (NRF). 5* is set to zero, and we identify the effective potential V with the zeroth component Uy, which 
leads to 



H NRF = ^/m 2 + P nrf + V - ( 43 ) 
The scalar potential U in any frame is now defined as 



u = \I H nrf - Pnrf - m - ( 44 ) 

The effective potential V , which describes many-body interactions of the nucleons, can be parametrized as a sum of 
a Skyrme part, only depending on the density p, and a momentum dependent part. We take a parametrization of 
Ref.|i, 

P(r) , D MOV , 2C f dV f(r,p) 



Po \ Po J Po J (27r) 



V(f,p)=A^+B[^\ + —g * (45) 



f 



( p-p' V 
I A J 
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Throughout this work we use A = -29.3 MeV, B = 57.2 MeV, C = -63.5 MeV, r = 1.76 and A = 2.13 far 
These parameters are fitted to the saturation density of nuclear matter, and also to the momentum dependence of 
the nucleon optical potential as measured in pA collisions 68] . 

Note that the introduced quantities describe the structure of the nucleus and do not depend on the particular 
nuclear reaction. Therefore, there are no free nuclear parameters for vA scattering. 

If a particle is bound in the nucleus it acquires an effective mass, defined as 

M eff = M + U(f,p). (46) 

As consequence of the density and momentum dependence of the scalar potential the effective masses of initial and 
final particles are different even if their masses at rest are equal. This fact is taken into account for QE scattering by 
replacing M and M' in the cross section formula (Eq. (0) and in the QE hadronic tensor (Eq. JUJ) with 

M -> M N + U{r,p), (47) 
M' -> M N + U(f,p'). (48) 

Analogously, for A excitation, the effect of the binding potential is considered by substituting M and M' in the 
cross section formula (Eq. cf. Eq. (|30|l ) and in the A hadronic tensor (Eq. (|19(l ) by 

M -» M N + U(r,p), (49) 
M' -» M A + U A (r,p'). (50) 

We consider that the A is less bound in the nucleus than the nucleons by setting the A potential Ua to 2/3 of the 
nucleon potential U. This choice is motivated by the phenomenological value of Ua — —30 MeV at normal nuclear 
density po [35| . 

Furthermore, we take into account that the width of the A is modified in the nuclear medium. While the nucleons 
inside the nucleus are constrained to have momenta below the Fermi momentum, there is no such constraint for the 
production of the resonances. Their decay, however, is influenced by Pauli blocking, e. g. a resonance decaying into 
a pion nucleon pair is Pauli blocked if the nucleon's momentum is below the Fermi momentum. Therefore, the width 
of the resonance inside the nuclear medium is lowered due to Pauli blocking. On the other hand, it is increased by 
collisions inside the medium. The collisional width T co u accounts for additional decay channels of the A inside the 
nucleus. Through two-body and three-body absorption processes like AA^ — > AW or ANN — * NNN, the A can 
disappear without producing a pion, while via AN — ► itNN additional pions can be produced. Also elastic scattering 
AJV — > AA^ contributes to the collisional broadening. For T co a we use the results of Oset and Salcedo 69] that have 
been conveniently parametrized. Therefore, we have for the total in-medium width 

r-f rf = f + r coH) (51) 

which then replaces V in Eq. I|30|l . Note that even though the Pauli blocked decay width T reduces the total in-medium 
width, the overall effect is a broadening of the A in the medium (cf. e. g. Ref. |7(|). Numerical results, which illustrate 
the influence on the cross section of the different in-medium modifications in comparison to the elementary case, are 
presented in Sec. IIII Al 

Notice that we do not yet include any modifications of the form factors due to the nuclear medium. Moreover, it 
has been shown for QE scattering by Singh et al. 0] and Nieves et al. 0, |2^| that the renormalization of the weak 
transition strength in the nuclear medium has a non-negligible impact on the cross section. We shall consider this 
additional modification in the future. 
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C. Final-state interactions 



Final-state interactions (FSI) of the produced particles — in our case nucleons and A resonances — are implemented 
by means of the semi-classical Boltzmann-Uehling-Uhlenbeck (BUU) transport model. This appro ach allows a full 
coupled-channel treatment of the FSI. Originally developed to describe heavy- ion collisions |72l FH. l74| . the Giessen 
BUU model (GiBUU) has been extended to describe interactions of pions, photons and electrons with nuclei 1^3, 
I73.l73.l77| . In particular, it has been successfully tested against the existing data on the interaction of photons and 
electrons with nuclei |78i l79j | — an important prerequisite for any model aiming at the description of the interaction 
of neutrinos with nuclei. The applicability to and test in many different nuclear reactions is a strength of our model, 
which describes many different processes using the same physics input. In this respect, the study of vA scattering is 
a natural extension of our previous work and does not introduce any new free parameter. Here we review only the 
main ingredients; details of the Giessen BUU model can be found in the references given above. 

The BUU equation describes the space-time evolution of a many-particle system under the influence of a mean-field 
potential and a collision term. A separate transport equation for each particle species is required. For a particle of 
type i it is given by 



[d t + V p H • V r - V r H • VpJ Fi(f, p,m;t) = Icoii[F h F N ,F„,F A , ...], 

where Fi( r , p, m; t) is the phase-space density at time t; f , p are the coordinates and the momentum of this particle. 
H is the relativistic Hamilton function of a particle of mass m in a scalar potential U, which is given by 



H= \J[m + U(r, p)] 2 + p 2 . (52) 



For nucleons, the scalar potential U was explicitly discussed in Sec. IIIBI The same potential is used for all baryons 
except the A resonance, for which we take Ua = \V ■ We include neither potentials for the mesons nor Coulomb 
corrections. 

Icoii is the collision integral which accounts for changes in the phase-space density due to elastic and inelastic 
collisions between the particles. If the respective particle decays into other hadrons, this also accounts for changes in 
Fi. Thus, I con consists of a gain and a loss term |68l l77| which include also Pauli blocking. All BUU equations are 
coupled through I co u and, with less strength, also through the potential in H. In this coupled-channel treatment of 
the FSI, our model differs from standard Glauber approaches since the collision term allows not only for absorption 
but also for side-feeding and rescattering. 

In between the collisions, all particles are propagated in their mean-field potential according to the BUU equation. 
Our model allows for off-shell transport of resonances as described in Refs. |68H77j. 

In the non-strange mesonic sector we include tt, r\, p, lo, a, <p, 77'. A full list of all mesons (also those containing 
strangeness and charm) with their properties is given in Ref. |77|. Besides the nucleon and the A resonance, we include 
in the baryonic sector 29 other nucleon resonances |6Sj . Their properties are taken from an analysis done in Ref. |63| . 
They can couple to the channels Nit, Nrj, AK, Nu, Att, Np, No, 7V*(1440)7r and Ap. Also in the baryonic sector, 
strange and charmed hadrons are included and we refer to Refs. |68l l77j for a full list. Finally, we take into account 
the corresponding antiparticle of each baryon. More details on the cross sections in the resonance region and their 
in-medium modification are given in Refs. [HUH El- Above invariant energies of y/s > 2.6 GeV for baryon-baryon 
collisions or y/s > 2.2 GeV for meson-baryon collisions, the FRITIOF string fragmentation model is used [8flj . Thus, 
our transport model is also able to describe the final-state interactions of high-energy particles and can therefore be 
extended in the future to the study of neutrino scattering in the DIS regime, where it has already been successfully 
applied to the electroproduction of hadrons in nuclei [TtJ ■ 

However, the most important particles for this work are the nucleon, the A resonance and the pion. For the NN 
cross section and its angular dependence we use a parametrization of Ref. [8l| which has been fitted to available data. 
The A resonance is propagated off-shell in our model — its decay is isotropic. Absorption processes of the A are 
implemented by an absorption probability depending on the in- medium width given in Eq. I|51|) . The treatment of 
pions, and the whole ttNA dynamics, in our BUU model has undergone numerous previous tests in A A— > ttX |72j . 
tt A — > X [Tol 182^ and ^ A — > tt X [Tg reactions (see Ref. [jjj for the cross sections used in the present article). In 
particular, in Ref. [8j| quantitative comparisons of calculated 7 A — > tt° X cross sections to experiment have been 
given. These results, which contain the side-feeding of the it channel from originally produced charged pion channels, 
show that the measured pion momentum distributions are described very well by our theory. The method has, because 
of its semiclassical nature, a lower limit of validity of about 20 - 30 MeV pion kinetic energy [jjj; this limitation is 
of no concern for the results to be discussed later. The comparison in Ref. j^l to data obtained with a few hundred 
MeV photon beam, roughly corresponding to the energy regime treated in this article, also shows that the deviations 
from experiment are typically of the order of 20%. We thus expect a similar systematic uncertainty in the pion results 
reported in this article. 
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As we shall see in Sec. lIIIBl FSI lead to absorption, charge exchange, a redistribution of energy and momentum as 
well as to the production of new particles. 



III. RESULTS 



We now present different results for CC neutrino-induced reactions on nuclei at neutrino energies up to 2 GeV 
where we have studied, as a representative example, the reaction 56 Fe — > fx~X. Iron will be one of the targets of 
the upcoming experiment MINERVA [f| which aims at a precision measurement of the neutrino cross sections. In 
order to compare with other calculations, some results are also presented for scattering off 16 O and 40 Ar. However, 
we want to emphasize that our model is applicable for CC and NC scattering of v and v of all flavors off any nucleus 
(from 12 C on). 

In this section we first discuss the inclusive cross section for QE scattering and A production and its sensitivity to 
the in-medium modifications introduced in Sec. IIIBI By accounting for FSI within the Giessen BUU model we can 
investigate in detail pion production and nucleon knockout — results are presented in Sec. IIII Bl 



A. Inclusive cross sections 



By inclusive we understand the vA reaction in which a muon is produced, regardless of the rest of the outcome. 
The cross section for the vA reaction is given by an integration over the contributions of all nucleons. 

In Fig. 01 we show the inclusive double differential cross section on 56 Fe at E v = 1 GeV taking into account 
Fermi motion, Pauli blocking, binding potentials and in-medium widths. At Q 2 < 0.4 GeV two peaks can be clearly 
distinguished. The one at higher corresponds to quasielastic events, whereas the one at lower results from A 
production. At higher Q 2 , the two peaks overlap, the distinct peak structure vanishes and the inclusive cross section 
tends to zero. 

To understand the influence of the different in-medium effects, we cut the double differential cross section at a fixed 
value of Q 2 = 0.15 GeV 2 (see Fig. EJ). The solid line corresponds to the full model calculation as in Fig.0J 

First, we want to draw attention to the QE peak (right peak in Fig.EJ. The line denoted with "elementary" is the 
vacuum cross section which can be obtained from Eq. In the case of QE scattering, the differential cross section 
is a (5-function in E^ whose position is marked by an arrow in Fig. [SJ This position, at which QE scattering occurs, 
is given by 



_Q 2 
2M N ' 



E » = E »- T±F-_- (53) 



If we "switch on" Fermi motion, denoted as "+ Fermi" in Fig. [S] (it overlaps with the dash-dotted line), we find that 
QE scattering is now possible for a range of values of with Q 2 and E v fixed. This can be understood from the 
on-shell condition 

s = (k + pf = Ml -Q 2 + 2E q E - 2p ■ q = M 2 N , (54) 

which has multiple solutions for E^ — E v — E q due to the Fermi momentum p of the nucleons. Pauli blocking 
(dash-dotted line in Fig. O causes a reduction of the cross section when E^ gets close to E v . The effect is actually 
more apparent when the neutrino energy is smaller. Furthermore, we take into account the binding of the nucleons in 
a mean-field potential (short-dashed line, at the QE peak this line overlaps with the solid line) by introducing effective 
masses as outlined in Sec. Ill Bl The inclusion of the potential lowers and shifts the peak to smaller energies. 

Turning to the A excitation (left peak in Fig. |SJ), we obtain the vacuum cross section A x d 2 a V N /(dQ 2 dE ll ) from 
Eq. J2J (denoted by the long-dashed line) by replacing the J-function according to Eq. I|30|) which then leads to the 
peak structure at lower E^. For the same reason as for QE scattering, we obtain a broadening of the peak when 
Fermi motion is included (dotted line, which coincides with the dash-dotted line). The inclusion of Pauli blocking 
has of course no impact on the A production cross section. Taking into account the nuclear binding of the initial 
nucleon and also the potential of the final A, we observe a broadening and a shift of the peak. Finally, we include 
the in-medium modification of the width of the A resonance, labeled in Fig. [S] as "+ in-medium width". The vacuum 
width is replaced by a sum of the vacuum width, modified due to Pauli blocking, and a collisional width accounting 
for additional channels in the medium, as discussed before. Since the cross section scales with the inverse of the width 
this also lowers the peak. 

The integration over leptonic degrees of freedom yields the differential cross section da/dQ 2 and the total cross 
section a. They are shown in Fig.|S|for QE scattering and in Fig.Qfor A excitation. In both figures, the in-medium 
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FIG. 4: Inclusive double differential cross section d 2 cr/(dQ 2 d-E Al ) on 56 Fe at E v — 1 GeV. The higher peak is due to QE 
scattering, the lower one is due to A excitation. Included are all in-medium modifications: Fermi motion, Pauli blocking, 
binding energies and the modification of the A width in the medium. The 3-dimensional plot of the upper panel is projected 
to the Q 2 — plane in the lower panel where the contour lines increase from to 680 x 10 _38 cm 2 /GeV 3 equidistantly by 
20 x 10" 38 cm 2 /GeV 3 per contour. 



modifications are included subsequently as indicated in the plots — their explanation was given above. In both cases 
the total cross section a is reduced in the medium (solid line) compared to the one for scattering on a free nucleon 
{(Jn] long-dashed line): a < A x er/y 

Our result for the QE total cross section can be compared with previous calculations on 56 Fe. Paschos et al. [84| 
only include a "Pauli" factor and neglect other nuclear effects such as Fermi motion and binding. In our calculation, 
the latter lead to a significant reduction of the free cross section that can not be obtained with Pauli blocking alone. 
Athar et al. [7l| consider, as we do, all these effects and further account for renormalization of the stren gth in the 
medium. Thus, they obtain a cross section even lower than ours (approx. 30% versus 10%). Juszczak et al. .85] apply 
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FIG. 5: Inclusive double differential cross section d 2 a/(dQ 2 dS M ) on 56 Fe at E v = 1 GeV and Q 2 = 0.15 GeV 2 . The left peak is 
due to A excitation, the right one to QE scattering. The dashed line shows the elementary cross section for A production. The 
position of the <5-function of the QE cross section is indicated by the arrow. We subsequently include Fermi motion (dotted 
line, this line overlaps with the dash-dotted line except at E^ — > 1 GeV), Pauli blocking (dash-dotted line), binding energies 
(short-dashed line) and the in-medium modification of the A width (solid line). 




Q 2 [GeV 2 ] E v [GeV] 



FIG. 6: Inclusive QE cross section on 56 Fe. The left panel shows the differential cross section da/dQ 2 at E v = 1 GeV, the right 
panel the total cross section a as a function of E v . The in-medium modifications as labeled in the plot are added subsequently 
to the calculation. The solid line shows the full model. "Elementary" denotes the vacuum cross section multiplied by the mass 
number A. 

a Fermi gas in the local density approximation with Pauli blocking. Furthermore, they use a momentum dependent 
binding potential. They find — as we do — a reduction of the free cross section due to the Fermi momenta and Pauli 
blocking of the nucleons and a further decrease due to the binding. 

The total cross section for A excitation can be compared to Athar et al. [ZlJ who find a reduction of 5 - 10% with 
respect to the elementary ones. As in our model the authors have included in-medium modifications of the width of 
the A but not any nuclear binding. However, through the consistent inclusion of binding energies, we obtain an even 
bigger reduction of about 20%. 

Finally, we show in Fig. [H] the differential cross section da/dQ 2 at a neutrino energy of 1 GeV (left panel) and 
the total inclusive cross section a as a function of E v (right panel). The solid line represents the full model with all 
in-medium effects as described above. We plot also the contributions from QE scattering and A excitation which add 
up to the full result. 

Apart from H2 and D2 scattering data from ANL and BNL discussed in the previous section, data were also 
taken for quasielastic and inelastic scattering on heavy (Neon, Propane and Freon) targets e. g. at Gargamelle at 
CERN [861 |87j and at the Serpukhov bubble chamber SKAT |88|, |89( . These experiments were limited in statistics 
with large neutrino flux uncertainties. Unfortunately, the large error bars do not allow to draw any conclusion about 
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FIG. 7: Same as shown in Fig. [U for A excitation. 
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FIG. 8: Inclusive cross sections on Fe. 



the role of nuclear effects on those measurements. For this reason we do not show them here. In any case, our 
results for the total QE cross section agree with these data within their error bars. In the case of pion production on 
heavy nuclei, there are scarce data points only at energies above the region considered here. The data will become 
more sensitive to in-medium processes with the cross section experiments MINERVA and FINeSSE [jj with such 
a precision measurement of both the quasielastic and the inelastic neutrino-nucleus cross sections, that a meaningfull 
comparison with theoretical models will be possible. 



B. Exclusive channels 



In the case of exclusive channels, even if the cross section for the vA reaction is again given by an integration over 
all nucleons, for each distinct channel, e. g. 7r+ production, the contributions have to be weighted with the multiplicity 
of the given final state. This multiplicity is calculated with the coupled-channel BUU transport model. 



1. Pion Production 



We start our discussion of pion production with the double differential cross section shown in Fig. and Fig. EH 
There we plot the exclusive cross section for neutrino-induced ir + and ir° production on 56 Fe at E v = 1 GeV as a 
function of and Q 2 . In this calculation we included all medium modifications of the elementary cross section as 
explained in Sec. Ill Bl The left panels shows the results without FSI while in the right panels they are included. 
Once produced the A can decay or interact via A7V -> NN, ANN -> NNN, AN ->■ irNN or AiV -> AN. The 
produced pions interact through irN — > irN, NNw — > NN and irN — > A, i. e. they can scatter elastically, undergo 
charge exchange or be absorbed. This results in the creation of additional pions or their absorption. 

The cross sections in Fig.[S]and Fig.^|peak clearly at the position of the A (compare to the lower panel of Fig.^J 
which indicates that most of the produced pions come from the initially produced A. Pions can also originate in 
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FIG. 9: Double differential cross section da/(dQ 2 d£' fl ) for ty + production on 56 Fe at E v = 1 GeV. The cross section is mapped 
to the Q 2 — Efj, plane. In this result all in-medium modifications of the elementary cross section are included. The right 
panel additionally includes FSI. The contour lines are equidistant every 5 x 10 _38 cm 2 /GeV 3 from to 95 (left panel) and 
35 x 10~ 38 cm 2 /GeV 3 (right panel), respectively. 
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FIG. 10: Same as shown in Fig.EJfor 7r° production. The contour lines are equidistant every 1 x 10 38 cm 2 /GeV 3 from to 21 
(left panel) and 10 x 10 _38 cm 2 /GeV 3 (right panel), respectively. 



multistep processes initiated by QE scattering. There, the produced nucleon can scatter in the nucleus and create 
pions through NN — > NA or NN — > NNtt. However, this effect, which can occur only at high momentum transfer 
Q 2 , is less important as can be seen in Fig. El There, we show the total cross section for tt + (left) and tt° production 
(right panel). The dashed lines show the result without final-state interactions, while the solid lines denote the result 
of the full calculation. Furthermore, we show the contribution from initial A excitation (dotted line) and from initial 
QE events (dash-dotted). The cross sections without FSI for ir° production is significantly lower than the one for ir + . 
This difference is a consequence of the primary interaction mechanism: 

vp -> rA++, (55) 
vn -> rA+. (56) 
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FIG. 11: Total cross section for 7r + (left) and ir° (right panel) production on 56 Fe. The dashed line shows the results without 
FSI (only the decay of the A is possible); the results denoted by the solid line include FSI. Furthermore, the origin of the pions 
is indicated (QE or A excitation). 



The first process is enhanced by an isospin factor of three. These As decay into pions by 

A++ -> P7T+, (57) 
A+ -> pn", nn + . (58) 

i r 2 i — 1 

With the isospin amplitudes of these processes, we obtain a ratio of tt + : n = Z + (l/3) 2 iV (-V2/3) N\ = 

4.4 : 1 (JV and Z are the proton and neutron numbers) for the cross sections without final-state interactions. FSI, 
however, change this ratio. Indeed, the comparison of the ir + channel to the ir° channel (left and right panel of 
Fig. Ill|) reveals big differences. For tt + we find a strong reduction of the cross section due to FSI, while in the 7r° 
channel this reduction is much smaller (compare dotted and dashed lines). The additional strength in the ir° channel 
is a consequence of the "disappearance" in the dominant n + channel: 7r + undergo charge exchange reactions like 
TT + n — > ir°p contributing in this way to the ir° channel. This leads to the observed side-feeding. Side-feeding in the 
opposite direction is strongly suppressed by the ratio of ir + to 7r° production on the nucleon. 

As said before, ir° and tt + production through FSI of QE scattering is not very sizable and happens only if the 
neutrino energy is high enough. However, the effect is relatively more important in the ir° channel than in the ir + 
one. This follows from the fact that, while the production of both ir° and ir + from initial quasielastic scattering is 
basically the same, the ir° : ir + ratio from initially produced A resonances is roughly a factor of 4 smaller as just 
outlined. Thus, this effect also enhances the ir° channel due to FSI. 

7T _ cannot be produced directly in the neutrino-nucleon reactions, but only via final-state interactions. Thus, they 
play only a minor role as can be seen in Fig. 1121 where we plot the total cross section for tt + , tt° and ir~ production 
including FSI. Note that this situation is reversed in antineutrino reactions where only A~ and A can be produced 
in the initial interaction. 

Further details can be brought up by studying the pion kinetic energy distributions. They are shown in Fig. 1131 
(Fig. I14|) for 7r + (n ) production at different values of E u . The dashed lines show again the result without final-state 
interactions and the solid lines the result of the full calculation. The contributions from initial A excitation (dotted 
line) and from initial QE events (dash-dotted) are also plotted. ir° and tt + production through FSI of QE scattering 
contributes mostly to the low energy region of the pion spectra because of the energy redistribution in the collisions. 

The maximum of the solid curve (i. e. the calculation with final-state interactions) peaks at 0.05 - 0.1 GeV in all 
cases shown in Fig. 1131 and in Fig. [21 This is due to the energy dependence of the pion absorption. The absorption 
is higher in the resonance region where the pions are mainly absorbed through the reaction nN — > A, followed by 
AN — > NN. This strong reduction for high energy pions and the corresponding shift of the maximum to lower 
energies can be seen by comparing the dashed and the solid lines. These absorption processes equally affect 7r + and 
7T° yields. But pions do not only undergo absorption when propagating through the nucleus. Of particular importance 
for pions of all energies is elastic scattering ttN — > irN which redistributes the kinetic energies, again shifting the 
distribution to lower energies. It is important to stress that similar patterns are obtained within our BUU model for 
7T° photoproduction in nuclei in a good agreement with data as can be seen in Fig. 14 of Ref. [83 ■ 

The different scale of tt + and 7r° (Fig. ^| and Fig. I14|) is a consequence of their different production rates in the 
neutrino-nucleon reaction. This leads to the already discussed side-feeding from the dominant tt + channel to the 7r° 
channel. Also, pions produced from initial QE events, contribute relatively more to the ir° channel. For this reason, at 
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FIG. 12: Total cross section for ir + (dashed), tv° (solid) and 7v~ (dash-dotted) production on 56 Fe versus the neutrino energy. 
These results were obtained with a full-model calculation. 

E v > 1 GeV we even get an enhancement of ir° at low kinetic energies compared to the calculation without final-state 
interactions. 

In Fig. 1151 and Fig. 1161 we show the angular distribution for 7r + and 7r° production. Notice that the forward 
peaked tt + are absorbed and thus a flatter distribution is obtained. Pions produced in FSI from initial QE events do 
not have any preferred direction, which leads to a rather flat contribution (dash-dotted line). The tt° channel gains 
strength from the explained side-feeding, therefore, the absorption is less and the cross section is even enhanced for 
large scattering angles. 

Finally, we compare with calculations available in the literature. In Ref. I^j, the FLUKA cascade model predicts 
a large pion absorption in the A region. Namely, for 1 GeV energy they find that only 55% of the charged pions 
leave a Fe nucleus. We obtain a very similar result of 51% as can be seen in the left panel of Fig. El Large in-medium 
effects are also reported by Singh et al. [2j| . They obtain an overall reduction of 40% in the CC A production of 
pions at E v — 0.75 GeV on 16 0. By comparing our free cross section shown in the right panel of Fig. [7] with the sum 
of all pion cross sections of Fig. ^] at the same energy we find a reduction of 67% on a considerabl y h eavier nucleus 
( 56 Fe). For 16 we obtain approximately 50%. Using the model of pion rescattering of Refs. [30l l34| Yu [9l| finds 
for 16 O a total reduction of 40 - 60% in the ir + channel and - 25% in the ir° channel depending on the absorption 
model. The smaller reduction in the tt° channel is also due to side-feeding from the dominant tt + channel followed 
by charge exchange. At 2 GeV, we obtain an overall reduction on 56 Fe of 43% and 2% for ir + and 7r°, respectively 
(compare dashed and solid line of Fig. Illfl . In the case of 16 we find also at 2 GeV a reduction of 28% for ir + and 
an enhancement of 11% in the ir° cross section. 

Here we shall recall that all the results obtained for pion production carry an uncertainty related to the poor 
knowledge of the A form factors as discussed in Sec. Ill A"2*l It is also clear from Fig. [31 that the contribution from the 
higher resonances will affect the pion yields. In the tt + channel we expect a minor impact of the order of 10% since 
most of them come from the A dominated ir + p channel, while there could be a change of up to 30% in the ir° channel 
at higher energies due to the smaller cross section. We shall account for this missing strength in the future. 

2. Nucleon knockout 

The channel under investigation now is nucleon knockout. We take into account all nucleons which leave the nucleus 
due to the vA reaction. In Fig. 1171 (Fig. I18fl we plot the exclusive cross section for proton (neutron) knockout on 56 Fe 
as a function of Q 2 and for E v — 1 GeV. All medium modifications as explained in Sec. Ill Bl are included. In 
the right panels we additionally take into account all FSI. Two clearly separated peaks can be seen at low Q 2 in both 
figures. The one at higher E^ is due to the nucleons produced in initial quasielastic events, whereas the one at lower 
E^ results from A production. This clear separation is lost at higher momentum transfer. We have seen before that 
the inclusive cross section is smeared out with increasing Q 2 due to Fermi motion; the shape of the inclusive cross 
section (Fig. is reflected in the exclusive proton and neutron production cross section plotted here. 

The significant difference in scale for proton and neutron knockout is generated by the initial neutrino-nucleon 
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FIG. 13: Kinetic energy differential cross section for ty + production on 56 Fe versus the pion kinetic energy T n at different values 
of E v integrated over Q 2 and The dashed line shows the results without FSI interactions (only the decay of resonances is 
possible), the results denoted by the solid line include FSI. Furthermore, the origin of the tv + is indicated (QE or A excitation). 



production process: QE scattering produces only protons via 

vn — > l~p; (59) 

no neutrons are produced in the initial quasielastic interactions (cf. left panel of Fig. I18fl . Also the A production 
mechanism favors protons 

vp^> /"A++ -> rpn + , (60) 

vn-> Z"A + -> rpn°, rnvr + , (61) 

since the first process is enhanced by a factor of three. Using isospin amplitudes we obtain a ratio of p : n = 



— a/2/3) N (1/3) 2 N — 9.8 : 1 (N and Z are the proton and neutron numbers) in the A region for 56 Fe. 

Therefore, in a calculation without final-state interactions, proton and neutron knockout differ by about a factor of 
ten in the A region. 

With final-state interactions, this scenario changes. On the right panel in Fig. ^] we can see that neutrons are 
produced in the quasielastic peak region through multistep processes. The initial neutrino-nucleon QE reaction 
produces only protons, which undergo elastic and inelastic nucleon-nucleon collisions in the medium via NN — > TV TV, 
NN — > NNir or NN — > TV A. This leads to charge exchange, and thus, to neutron production. The situation is 
similar in the region of the A. 

The total cross sections for proton and neutron knockout are shown in Fig. 1191 The solid line, showing the result 
with all final-state interactions included, lies well above the one without FSI (dashed line) already for the protons, but 
even more so for the neutrons. This enhancement is entirely due to secondary interactions and cannot be obtained in 
a Glauber treatment. Furthermore, it is indicated in Fig. 1 191 whether the knocked out nucleon stems from initial QE 
scattering or A excitation. In contrast to the pion case, both contribute to the total cross section, even though with 
different weight depending on the neutrino energy. This weight follows simply from the total inclusive cross section 
(cf. Fig. El) an d thus explains the small contribution of the A at E v — 0.5 GeV, which gets larger at higher energies. 
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In Fig.l2Uland Fig.[^we present the kinetic energy differential cross section for proton and neutron knockout versus 
the kinetic energy for different values of E v . The line styles are as in the previous figures. FSI strongly modifies the 
shape of the distribution. High energy protons rescatter in the medium. As a consequence the flux at higher energies 
is reduced and a large number of secondary protons at lower energies appear (cf. Fig. I2U|) . Also low energy neutrons 
are produced through charge changing FSI as can be seen in Fig. I2T1 where, in the case without FSI, the cross section 
almost vanishes. 

The angular distribution is plotted in Fig. 1221 and Fig. [231 the nne styles are as in the previous plots. In the case 
of the protons, it is clearly seen how FSI shift the strength from small (~ 30°) to large angles (> 60°) as expected. 
The angular distribution for the neutrons shows an overall increase as a result of the described side-feeding. 

Finally, we want to compare our results with others. Fig. 1241 shows our results together with those obtained by 
Nieves et al. [28|. Since in this calculation only QE scattering is included, we have "turned off" the A excitation. 
The difference between the results without FSI (dashed versus dash-dotted line) is due to the RPA correlations taken 
into account by Nieves et al. which cause a reduction of the cross section and a spreading of the spectrum allowing, 
for larger energy transfers. However, when the rescattering of the outgoing nucleons inside the nucleus is considered, 
both calculations lead to very similar cross sections, namely a reduction of the flux for higher energetic protons and 
to a large number of secondary low energy protons (solid versus dotted line). 

The flux reduction can also be achieved with optical potential models. However, by simple absorption of the 
nucleons they do not account for the rescattering in the medium which leads to the large number of secondary 
nucleons. Nucleons are not just absorbed but — through rescattering — ejected with a different energy, angle and/or 
charge. In particular, with a model not accounting for FSI in a realistic way, there would be no neutrons in the final 
state when the initial vN collision was quasielastic. 



IV. SUMMARY AND CONCLUSIONS 



In this work we have investigated neutrino interactions with nuclei. The model presented here is able to describe 
neutrino reactions on both nucleons and nuclei at intermediate energies on the same basis as photonuclear quasielastic 
and meson production processes. 
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FIG. 15: Angular distribution for 7r + production on 56 Fe at various neutrino energies. The angle 9 is measured with respect to 
the direction of the incoming neutrino. Labels are as in Fig. 1131 and Fig. 1141 



We have focused on the region of the quasielastic and A(1232) peaks. For vN collisions a fully relativistic formalism 
is used. The extension to finite nuclei has been done in the framework of a coupled-channel BUU transport theory 
where we have taken into account in-medium modifications due to Fermi motion, Pauli blocking, nuclear binding and 
collisional broadening of resonances. This gives us the possibility to study exclusive channels taking into account 
in- medium effects and final-state interactions. 

To summarize our results on neutrino-induced pion production, we have found that the in-medium effects and, 
especially, final-state interactions reduce the exclusive pion cross section and also give rise to a small fraction of 7r~. 
Quasielastic scattering followed by it production in NN collisions contributes only weakly to the pion production 
cross section. In the kinematical region under investigation, the pions originate mainly from the initial A excitation. 
Furthermore, we have found an enhancement of the ir° channel through side-feeding from the dominant ir + channel. 

For nucleon knockout, we have found that the influence of the final-state interactions is significant. High energy 
nucleons rescatter in the nucleus, which leads to a decrease of the flux at higher energies, but also to a large number of 
secondary nucleons at low nucleon energies. While the elementary quasielastic reaction cannot produce neutrons, but 
only protons, we found that, as a consequence of the final-state interactions, a large fraction of neutrons is nevertheless 
produced. Also in the A region we found a large enhancement of neutrons due to final-state interactions. 

We conclude that in-medium effects in vA scattering, and in particular FSI, are essential. Their understanding — 
within a well tested transport model — is crucial for current and future experiments. 

At the point of submission, we became aware of a article by Cassing et al. |92l | studying the problem of neutrino- 
induced pion production. 
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FIG. 16: Same as shown in Fig. 1151 for it production. 
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FIG. 23: Same as shown in Fig. 1221 for neutron knockout. 




FIG. 24: Muon neutrino-induced proton knockout on 40 Ar at E„ — 0.5 GeV by QE scattering. The dotted and the dash-dotted 
line show the results of Ref. |2S|, the solid and the dashed line are the results of our model without A excitation. 



